Interaction effects on quasiparticle localization in dirty superconductors 
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We study how quasiparticle interactions affect their localization properties in dirty superconduc- 
tors with broken time reversal symmetry - for example in a magnetic field. For SU(2) spin-rotation 
invariant (class C) systems, the only important coupling is the spin-spin triplet interaction, which 
we study within a renormalization group approach. Either an additional Zeeman coupling or a com- 
plete breaking of spin rotation symmetry renders all interactions irrelevant. These two situations 
realize, respectively, the non-interacting unitary Anderson and the "thermal" (class D) universality 
class. Our results imply a stable metallic phase in 2D for class D. Experimental implications are 
discussed. 
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The interplay between disorder and interactions in 
electronic systems underlies several interesting phenom- 
ena in solids. A famous, though poorly understood, 
example is the metal insulator transition in three- 
dimensional dirty solids. Although much theoretical 
progress has been achieved in understanding the Ander- 
son localization transition of non-interacting electrons, 
direct contact with experiments has been problematic 
due to the complicating effect of interactions. Interest- 
ing field theoretic attempts have been made to incorpo- 
rate interactions into the successful scaling theory of the 
Anderson transition M. However the resulting theories 
are complicated to analyse, and in most cases have thus 
far not led to a good description of real metal-insulator 
transitions. In two dimensional systems (2D), the sit- 
uation is even worse: even the possible stability of a 
genuine metallic phase in the presence of disorder and 
interactions is a matter of considerable current debate 
H . Recent work J3|,|| has made some progress in under- 
standing the field theory of the disordered, interacting, 
two-dimensional electron gas. 

As recently emphasized, the dynamics of quasiparti- 
cles in a superconductor provides a new, though still ex- 
perimentally relevant, context to address localization is- 
sues |p|-|l0||. All disordered superconductors fall into one 
of two categories according to the nature of their quasi- 
particle transport properties - superconducting "thermal 
metals" (with delocalized quasiparticles) or "thermal in- 
sulators" (with localized quasiparticles). Interesting dif- 
ferences arise with localization physics in a normal metal 
due to the lack of conservation of the quasiparticle elec- 
tric charge in the superconductor. 

In this paper, we take up the task of describing quasi- 
particle localization inside superconductors in the pres- 
ence of both disorder and interactions We argue 
that this problem is simpler than the corresponding prob- 
lem in normal metals. We explicitly identify physical 
situations in which the interaction effects are unimpor- 



tant for the long distance physics. The corresponding 
experimental systems thus provide a clear opportunity 
to study Anderson localization transitions, unhindered 
by interaction effects. We also demonstrate the stability 
of a thermal metal phase in 2D in the presence of both 
interactions and disorder inside superconductors under 
appropriate conditions. The insights gained from study- 
ing the superconductor may be valuable in developing an 
understanding of the normal metal. 

Within the standard mean field treatment of pair- 
ing, the dynamics of non-interacting BCS quasiparticles 
is governed by a quadratic Bogoliubov-deGennes (BdG) 
Hamiltonian, subject to static disorder in the normal and 
pairing potentials. The localization physics of quasipar- 
ticles has been studied previously in this approximation. 
A total of four universality classes, different from the 
three known standard classes for normal metals, have 
been found for the possible localization behavior of non- 
interacting quasiparticles inside dirty superconductors 
on length scales much larger than the mean free 
path. These correspond to BdG Hamiltonians with or 
without spin-rotation and/or time reversal symmetries. 
Since charge is not conserved in a superconductor, the 
nature of a phase (metal, insulator, or Hall insulator) 
manifests itself not in charge transport, but rather in 
thermal-transport or, when spin is conserved, in concomi- 
tant spin-transport. 

Here we investigate the effect of quasiparticle inter- 
actions, focusing on situations lacking time reversal in- 
variance. One simplication offered by the superconduc- 
tor (as compared to a normal metal) is that the long- 
range Coulomb interaction is always screened out by the 
condensate. Thus, the quasiparticle interaction is short- 
ranged. Moreover, lack of charge conservation renders 
the singlet density-density interaction unimportant alto- 
gether. As in normal metals, systems with broken time 
reversal symmetry offer the further simplification that in- 
teractions in the Cooper channel are also unimportant. 
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The most significant interaction then is a short-ranged 
interaction between the quasiparticle spin densities. 
Possible experimental realizations include certain 



heavy fermion superconductors 12 1, which typically have 
strong spin-orbit (S.O.) scattering, and superfluid He-3 
in porous media Q ; these fall into class D of || . 

Perhaps the most promising prospect is Type II super- 
conductors in strong magnetic fields [|| . It has been sug- 
gested that a "thermal insulator-metal" transition, which 
can be probed by ultralow temperature heat transport 
measurements, may be driven in such a superconductor 
by simply changing the magnetic field. If the Zeeman 
coupling to the magnetic field and S.O. scattering can 
be ignored, the transition is in the universality class de- 
scribing a superconductor with full SU(2) spin symme- 
try, but without time reversal symmetry (class C of Ref. 
^|). However, typical Zeeman energies are expected to 
be much bigger than the low temperatures necessary to 
clearly extract the electronic heat transport. Thus, the 
asymptotic critical properties will be described by a the- 
ory that includes the Zeeman coupling. In the absence of 
S.O. scattering this situation is formally in the same uni- 
versality class as that of spinless electrons in a magnetic 
field §. Short -range interactions are known to be irrel- 
evant in the latter model (lj] - a physical consequence 
of the Pauli principle. Therefore, the field-driven ther- 
mal metal-insulator transition in a bulk 3D type II su- 
perconductor provides an excellent (and possibly unique) 
opportunity to experimentally study the non-interacting 
3D (unitary) Anderson transition, in contrast to nor- 
mal metals where the long-range Coulomb interaction 
changes the universality class. Despite being irrelevant, 
the short-range interactions in this system can still affect 
the low frequency, finite temperature (T) dynamics: a 
calculation along the lines of Jl|| shows that the ther- 
mal conductance k for T — > behaves as k/T oc T 8 ( d ~ 2 \ 



Without interactions, 9 = the localization length ex- 
ponent. With interactions, 9 = min(v,^), where p is 
the "dephasing" exponent arising from the irrelevant in- 
teractions. To lowest order in the d = 2 + e expansion 
one finds p < 1v (p — 1.3 when naively extrapolated to 

d = 3.) 

If, on the other hand, even uniaxial spin rotation sym- 
metry is broken (class D of H) due to, for example, S.O 
scattering, then for the same physical reason as above, in- 
teractions are again expected to be irrelevant. In 3D, this 
is thus expected to realize the thermal metal-insulator 
transition discussed in ]7j]. In 2D, this implies that the 
(thermal) metallic phase in non-interacting models of 
such superconductors [|fP]lO|l is stable to the inclusion 
of interactions. This then provides a concrete theoreti- 
cal instance of a stable metallic phase in 2D, albeit for 
thermal transport. 

To study interaction effects in classes C and D, we con- 
struct non- linear sigma models [NLerM] which generalize 
those constructed by Finkclstcin for normal metals [Q. 
Within this framework we see the irrelevance of all inter- 
actions in class D, and that of singlet and Cooper inter- 
actions in class C. For the latter case we carry out a per- 
turbative Wilsonian renormalization group (RG) analy- 
sis: In 2 + e dimensions (e > 0) the remaining triplet 
interaction is found to be marginal, to 1-loop order, at 
the fixed point describing the (thermal) metal-insulator 
transition without interactions In 2D, we find that 
the triplet interaction strongly affects the weak localiza- 
tion correction to the spin- (and thermal-) conductivity 
in the metallic phase (at weak coupling), which changes 
sign for sufficiently attractive interactions. 

Class C: We start with a general BCS Hamiltonian 
for a dirty singlet superconductor, possessing spin rota- 
tion invariance, but no time reversal symmetry, in the 
absence of quasiparticle interactions: 
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J d d x ^ - M + V{x)j 4> a + (a(x)^\ + h.c.) 



(1) 



Here a =|, | labels the spin. The potential V(x), as well as the real and imaginary parts of the complex gap function 
A(x), are (static) random variables with Gaussian distributions of zero mean and (without loss of generality) of same 
variance. The vanishing mean of the complex A(x) is physically appropriate in the Type II superconductor in the 
mixed phase with randomly located vortices, and in other physical situations. We treat the disorder average using the 
Schwinger-Keldysh method (described for the unitary symmetry class Q - see also []l7| ) . Within this formalism, the 
fermions acquire an additional Keldysh index i = 1, 2 denoting the time-ordered and anti-time-ordered branches of 
the path integral. At zero temperature the Keldysh functional integral action corresponding to the above Hamiltonian 
is (spin and Keldysh labels supressed where possible) 



S = 



V 2 



ip(x,t) + irj / d d x— sgn(w)'(/' t (x, u)ip(x, u>) + 
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d d xdt ^{x,t)a, 

d d xdt V{x)^ (x, t)o- z ip(x, t) + A(x)V>} (a;, t)a z tp\ (x, t) + A* {x)ijj\ (x, t)a z i\)\ (x, t) 



(2) 
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where o~ z is a Pauli matrix in Keldysh space. The second 
term (77 > 0) specifies the time-ordering on the Keldysh 
contour. The symmetries are made explicit by introduc- 
ing a 4-component field on each Keldysh branch, 



Finkelstein form. In particular, the resulting Keldysh- 
Finkelstein action reads (after simple rescalings and uni- 
tary transformations): 



(x,t)\_J_( il>i(x,t) , 



Z = J [VQ]t 



S[Q]. 



S[Q]=S D [Q] + S int [Q] (4) 



«? rm- ( — 

satisfying the reality condition x\ = (CXi) T with C = ^d[Q\ — . 



(J, y Ty (Pauli matrices j2 a b act on the 'particle-hole' index 
a of Xi: an d ^afi on the spin-index.) When expressed 
in terms of the 8-component fermion field x(x,iX)\ the 
Keldysh action above is readily seen to be invariant under 
independent symplectic transformations U(lo) £ Sp(4) C 
11(4:) for each frequency u> (repeated indices summed), 



— TY(VQ) 2 + izTr{iua z - r?sign(u;))Q 
8ng 



Ty<J Z U 



when rj = 0. A finite rj breaks the symmetry down to a 
U(2) subgroup. By standard arguments the low energy 
degrees of freedom are described by a diffusion mode ma- 
trix field Q(x) = Q l jl a Jl (x) which decouples fermion bi- 

linears a % 2 (x iaa ( x i u{)Xjap(x, u> 2 )) , appearing after the 
disorder average. It carries spin (a,/3), Keldysh 
and frequency indices. Q takes values on the saddle point 
manifold specified in Eq. (Q) below, and its dynamics is 
governed by the Keldysh NLcrM action Eq. m). This 
formulation of the non-interacting disordered theory is 
entirely equivalent to those obtained using the replica 
trick or supersymmctry. 

Effect of interactions: In the absence of disorder, the 
interactions at the Fermi surface can be, most generally 
p8| , of density-density (singlet), spin-spin (triplet) and 
(singlet and triplet) Cooper types. Introducing three 
Hubbard-Stratonovich (H.S.) fields to decouple these 
four-fermion interactions, and condensing the H.S. field 
which decouples the Cooper interactions (gap function) 
in the spin-singlct channel, gives the BCS mean field the- 
ory. In the presence of disorder, but no interactions, this 
has the form of Eq. ([!]). Quasiparticle interactions arise 
from the (dynamical) fluctuations of the three H.S. fields. 
In particular, residual interactions in the Cooper chan- 
nel arise from the fluctuations of the amplitude and the 
phase of the gap function. The gapless phase fluctua- 
tions can be seen to decouple from the BdG quasipar- 
ticles at asymptotically low energies, and their only ef- 
fect is to render a long-range Coulomb interaction short- 
ranged This confirms the expected screening by 
the condensate, mentioned above. Averaging over disor- 
der in the presence of the three H.S. fields yields, along 
the lines of [|l6|, a theory of diffusion modes Q inter- 
acting with the latter. In the absence of fluctuations of 
the gap function amplitude, the remaining two H.S. fields 
can be integrated out straightforwardly. The singlet H.S. 
field is seen to decouple, and integration over the triplet 
H.S. field leads to an interaction term of the familiar 



Sint[Q] = iirUtz 



2 e< 



d d x 



(5) 



dtQffPQff (6) 



Here Tr denotes the trace over all indices of Q, includ- 
ing the frequency index u, and is the correspond- 
ing (double-) time Fourier transform. The saddle point 
manifold of massless modes of the NLcrM is described in 
frequency space by 

Q 2 = l, Q} = Q, Tr(Q)=0; t v Y, x Qt v 11 x = -Q T (7) 

where T, x exchanges positive and negative frequencies. 
Note that Qff transforms as a spin-triplet. Conse- 
quently, Eq.(^) is the only SU(2) spin-rotation invari- 
ant Finkelstein interaction term that can be written. 
This was expected: due to the lack of charge conser- 
vation there is no massless charge diffusion mode, and 
due to the lack of time-reversal symmetry, there is no 
massless Cooperon mode on the saddle point manifold. 
Hence Eq.(^|) represents the only non- vanishing interac- 
tion, which is that between the spin diffusion modes. 
Here 1/g is proportional to the (spin) conductivity, and 
Ut(< 0) is the (repulsive) triplet interaction strength. 

We perform a perturbative (Wilsonian) R.G. analysis, 
parametrizing the saddle point manifold by independent 
matrix fields V, similar to The wavefunction renor- 
malization is found by evaluating (Q), using the same 
logic as in the standard O(N) NLcrM The 1-loop 

R.G. equations, valid to lowest order in g (and z), but to 
all orders in UtZ, are found in d = 2 + e: 



dg = 
dl 2 
d(U t z) 
dl 
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2U t z 



log(l - 2U t z) 



g 2 + o(g*) 



dz 



0(g z ), -£= ~-{l + W t z)g z + Otf) 



(I is twice the log of the length scale). These are ex- 
tracted from the renormalization of the terms quadratic 
in V in the action. We also confirmed that the same 
R.G. equations are obtained from the renormalization of 
the terms of next higher order in V. This is a necessary 
condition for the renormalizability of the theory. 

Note that in contrast to other universality classes, Utz, 
a measure of the interaction strength, does not renormal- 
ize to 1-loop order in class C. 
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Without interactions (U t — 0), we recover the non- 
interacting fixed point at g = 5 + which is the ther- 
mal (or spin) metal- insulator transition considered in Q . 
We see from the second R.G. equation that the interac- 
tion U t z is marginal at this fixed point, at least to 1-loop 
order. Hence, showing that interactions are relevant or 
irrelevant would require working to higher-loop order. 

Next we specialize to 2D, focussing on the metallic 
phase where j < 1. Since Utz does not renormalize to 
lowest order in g , it will be constant over a wide range 
of length scales. Hence, in this range, the thermal- (k), 
and spin- (c s ) conductivity is given by the 1-loop result 

t 3h 2 K SV 



1 

Ait 2 



7-6 1- 



1 



217*2 



log(l - 2U t z) 



log 



Here CTq is the (bare) spin conductivity at the scale of 
the elastic mean- free path £ e , and a s its value at length 
scale L. This expression is perturbative in g = (^)l/a s , 
but UtZ need not be small. (Note that the quantity in 
square brackets goes to 1 as interactions are removed.) 
Thus, a repulsive (attractive) triplet interaction is seen 
to decrease (enhance) the weak localization correction to 
the (spin-, or thermal-) conductivity, as compared to the 
non-interacting case. For sufficiently strongly repulsive 
interactions, 2Utz < —0.37.., the g 2 -coefficient in dg/dl 
changes sign, implying that g flows back to weak cou- 
pling, at least until, potentially, higher-loop effects set in 
which may reverse the flow. 

Class D: Turning now to situations lacking spin- 
rotation as well as time reversal symmetry, we can see, 
without extensive calculations, that all Finkelstein-type 
interaction terms are absent. In this case, the diffusion 
mode matrix field Q l t 3 t (x) decoupling fermion bilinears 

a z (Xiaa( x > t i)Xjaa(x,t2)) carries only Keldysh and 
frequency, but no spin indices. Repeating the steps above 
Eq.(^), both the singlet and triplet H.S. fields are seen 
to decouple. This was expected since now, due to lack of 
spin-rotational symmetry, also the spin diffusion mode is 
absent from the saddle point manifold of massless modes. 
Indeed, the only possible Finkelstein interaction term as 
in Eq.(j^) (but without spin indices, a, (3) vanishes due 
to the antisymmetry of QJ* tAx), which follows from the 
reality condition below Eq.(||). 

Consequently, we conclude that the 3D thermal metal- 
insulator transition in class D j| is expected to be 
unmodified by quasiparticle interactions, as mentioned 
above. In 2D, in the absence of interactions, a (thermal) 
metallic phase, stable to quantum interference, can exist 
inside generic superconductors lacking spin rotation and 
time reversal invariances [Qj9|. (For spin-polarized ver- 
sions, see Ref. [|l0|.) The absence of any type of marginal 
Finkelstein interaction terms in class D, found above, 
thus implies that this 2D metallic phase remains stable 
even upon inclusion of any type of interactions. 



We end our discussion with a general caveat. The ques- 
tion as to whether or not additional potentially 'danger- 
ous', i.e. relevant or marginal, long-range (in time) ver- 
sions of the interaction terms are eventually generated 
upon the R.G., hinges upon a proper understanding of 
the renormalizability of the whole class of Finkelstein- 
type theories in general, which is currently lacking. How- 
ever, it is not expected that initially short-range interac- 
tions would generate such long-range ones, upon inte- 
gration over short-distance, short-time fluctuations in an 
R.G. transformation. These issues will be addressed in 
more detail in Ref. ]i~9| |. 

After our work was completed, and some of our results 
were reported in pd| ] , we learned about the work by Fab- 
rizio et al. |2l[], whose results for the correction to the 
2D conductivity and density of states for class C are in 
agreement with ours. 
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DMR-98- 76208 (C.C.). The work of T.S. was supported 
by the MRSEC program of the NSF through the grant 
DMR-98-08941. C.C. acknowledges the support from the 
Alfred P. Sloan Foundation. 
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